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Abstract 
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as 

5^ , We study the lensing of neutrinos by astrophysical objects. At the 

T^ I difference of photons, neutrinos can cross a stellar core; as a result 

the lens quality improves. While Uranians alone would benefit from 

Oh. this effect in the Sun, similar effects could be considered for binary 

r-| , systems. 
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1 Introduction 

In this note we want to investigate the possibihty of neutrino lensing by as- 
trophysical objects (stars, galaxies or rather galactic halo). At the difference 
of photons, neutrinos can cross even a star's core. As will be seen below 
this results in a much better focalisation. On the other hand, the price to 
pay is the extreme difficulty to detect neutrinos, and the comparatively poor 
angular resolution of "neutrino telescopes". For this reason, the only thing 
we can hope for is a signal intensification, rather than the spectacular pho- 
ton lensing patterns. The main equations are established in section 2, with 
details in the Appendix. Section 3 deals with neutrino interactions in a stel- 
lar medium, and specify which energy range can be studied in this way. In 
section 4 we consider a number of situations and establish the corresponding 
signal enhancement expected. Finally in section 5 we consider some practical 
examples. It appears clearly that the Earth - Sun distance is too small for 
a sizable effect to take place, (an observatory on Uranus would notice the 
enhancement of distant neutrino sources whenever they are aligned with the 
Sun). We then consider the case of galaxies, through their halo. We finally 
turn to binary systems. We will not in this paper comment on the possible 
origin of energetic neutrinos, but we remark that an interesting situation is 
met when one of the companions focuses the neutrino fiux originating from 
the other. 



2 Gravitational deflection of neutrinos 

In this section we study the defiection of neutrinos from straight-line motion 
as they pass through a gravitational field produced by a compact object 
of mass M. We will distinguish two cases: when the neutrino fiux passes 
far away from the object (OUTside solution), a situation equivalent to the 
gravitational lensing of photons, and, when the neutrino fiux passes through 
the object (INside solution). In the latter case, we consider three specific 
cases depending on the compact object density profile: constant density (an 
academic but constructive example), Gaussian distribution density (suitable 
for stars0), and Lorentzian distribution density (which could be associated 

^ The density profile of stars is not exactly Gaussian but we use it here in order to obtain 
simple analytical results. Such a description should be considered as a good approximation 



to a galactic haloQ). 

We begin by calculating the trajectory of a massless neutrino (or with 
mass very small compared to its energy) in the Schwartzschild metric under 
the assumption that M/r is everywhere small along the trajectory |jl|. The 
equation of the orbit is0 

d^ ^ 1 /^v 

* r-^li - * (1 - a?) ' 



r 



where h is defined as the impact parameter. Using the definition m = ^ 

dd) 1 



du Ji - m2 + 2Mm3 



(2) 



If we neglect the u^ term in Eq. (|^), all effects of M disappear, and the 
solution is 

rsin(0-0o) = & , (3) 

a straight line. In the limit Mu <C 1, or -Rsch. = 2M <^ 6, if we define 
y = u{l — Mu), Eq. (H) becomes 

-^^ l + 2My , ^,,^. 2 



Integrating Eq. (M gives 



+ 0{M'u') . (4) 



<PovT{y) = (j)o + — + arcsin(%) - 2mJ^ - y^ , (5) 

where the integration constant is defined as = 0o (with 0o the incoming 
direction) when the initial trajectory has r — »• cxd (or y — > 0). The particle 
reaches its smallest r when ^ = 0: 



dr „ L [^ 2M\ 62 ^ 1 . ^.,.2 2^ 



^ = i^Jl - (^1 - — j ^ = ^ 2/n.ax = ^ +0(MV) . (6) 



to the real case. 

"* The Lorentzian profile behaves as 1/r^ for large r, in agreement with velocity disper- 
sion curves for galaxies and clusters. 

^Geometrized units c = G = 1 are used throughout the paper. 



This occurs at the angle 

0out(i/ = 5) = '^o + ~ + 2 ■ ^^ 

It has thus passed through an angle ^ + | as it travels to its point of closest 
approach. By symmetry, it passes through a further angle of the same size 
as it moves outwards from its point of closest approach (see Ref. [Q). Then, 
the particle passed through a total angle of ^ + vr. If it were keeping to a 
straight trajectory, this angle would be vr, so the net deflection is 

. , 4M 

A0OUT = — • (8) 

For the case of a neutrino flux passing through the object, we must first, 
in order to study the neutrino trajectory, look for the form of the space-time 
in the region inside the object. Here, we restrict ourselves to the case of static 
spherically symmetric space-timesQ. In that case, the most general metric is 
(see Ref. [0] for details): 

where it is convenient to replace A(r) by 
1 



r 

For a static perfect fiuid|], Einstein equations imply 

^ = 47rrV(r) , 

(11) 

d<I>(r) m{r)+47rr'^p(r) 1 dp{r) 

dr r{r—2m(r)) p{r)+p{r) dr ' 



^Spherically symmetric space-times are reasonably simple, yet physically very impor- 
tant, since very many objects of importance in astrophysics appear to be nearly spherical. 
A static space-time is defined to be one in which we can find a time coordinate t with 
two properties: (i) all metric components are independent of t, and (ii) the geometry is 
unchanged by time reversal, t — > —t. 

^A perfect fluid in relativity is defined as a fluid that has no viscosity and no heat 
conduction in the momentarily comoving reference frame (MCRF) . A static fluid is a fluid 
that has no motion. 



where p(r), m{r) and p{r) are, respectively, the density, mass and pressure of 
the object at a radius r. For completeness, in the region outside the object 
we have p = p = 0, then 



m[r\ 



M = const. 



32${r) 



2M 



ds^ 



^]df 



+ 



dr^ 



-^ir 



+ r^dn"^ . 



(12) 



In the region inside the object, exact solutions to the relativistic equations 
are very hard to solve analytically for a given equation of state ^ . One inter- 
esting exact solution is the Schwarzschild constant-density interior solution, 
which we use here as an example of the framework needed to study other 
density profiles. 

Inside R, where R is the physical radius of the object, p 7^ 0,p 7^ 0. 
For a constant density profile p(r) = p, the equation of the orbit is (see the 
appendix for a detailed calculation) 



d(j) 
dr 
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(13) 



and the net deflection is 
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(14) 
where the outside solution is also included for completeness. 

Next, we analyze the solutions for the Gaussian and Lorentzian distri- 
bution densities. The Gaussian profile is a convenient approximation to the 
mass distribution in stars, while the Lorentzian profile is valid for galactic 
halos. In both cases, it is possible to neglect the pressure with respect to 
the mass density, p <^ p (see Ref. '^ for the so-called Newtonian stars), so 
we also have Anr^p <^ m. Moreover, the metric must be nearly flat, so in 
Eq. (p!OD we require m(r) <^ r. These inequalities simplify Eq. ( pT]) to 

d^{r) m{r) 



dr 



(15) 



and the metric in Eq. (^ to 

g2A(r) ^ ^ ^ 2n^ ^^^ g2$(r) ^ ^ ^ 2<l>(r) 

=^ ds^ = -(1 + 2<l>(r))dt2 + (i + Mr)) rf^2 ^ ^2^^2 



(16) 



For the case of a Gaussian density profile p{r) = p^e *" /'^o^ the equation 
of the orbit is (see appendix for details) 



dr 



M 



J.R^-r'^)lrl_^ 



X 



^ ro/-Re^^/''o^/2erf(R/ro)-l 
1 



'hi r2 



-I 1M rg/re ^/¥/2 erf (r/rp) ~1 

-R r2 /,.2 



(17) 



where the error function is defined as erf (2;) = ^ J^ dte * . The net deflec- 
tion is then 



4M 
b 



4M,i 



iib> R 



1 _ ''^ 



4A/ rp/iie^ /'•oV^/2 

^ ro/-Re^^'''"o0F/2erf(ii;/ro)-l 



(18) 



X 



1 - -^ erf (i?/ro) - e-^'/'-oeii Jl 



i?2 r-o 



if 6<i? 



It is very interesting (because we are close to a real case) to compare the pre- 
vious result with a naive approximation where for a given impact parameter 
b one studies the net deflection by a sphere of radius b and identical density 
profile 



A0 



approx- 



4m (b) 
b 



iM^(B?-r'^)/rl rp/re'' /'•oy¥/2erf(r/ro)-l 
-'^ ro/Re-^^/''o0F/2erf{_R/ro)-l 



for b < R 



(19) 



In Fig. |l|, we plot for comparison the exact result in Eq. ([18|) together with 
the naive approximation in Eq. (|l^). We have taken the Sun as example of 
a typical star. Then, the mass and the radius in Eqs. (|T8|,|l9|) are fixed to 
M = Mq = 1.48 km and R = Rq = 6.96 x 10^ km. The parameter tq is 
taken to be tq = 0.2Rq. As it is seen from Fig. |I], the maximal net deflection 



occurs at 6 ~ ro and is A0(6 ~ ro) — 3.2A0(6 = Rq), showing that the 
lensing effect inside the star is bigger than the outside effect (except for 
b < O.O4i?0). It is also shown that for the naive approximation the maximal 
deflection is displaced up to 6 ^ l-5ro and its value is around 20% smaller 
than for the exact result. So then, our calculation seems to be essential for 
a detailed analysis of the gravitational lensing of a neutrino beam by a star. 
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Figure 1: Net deflection A0 as a function of the impact parameter b for the 
case of a Gaussian density profile p{r) = poe~^ ^^o (see Eq. (|8])). The naive 
approximation of Eq. ( P^ ) is also included for comparison. In both curves, 
M = Mq = 1.48 km, R = Rq = 6.96 x 10^ km and ro = O.2i?0 are used. 



For the case of a Lorentzian density profile p{r) 






-TT^, the equation 
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(20) 



and the final defiection is 
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(21) 



iib<R 



Here, it is also interesting to compare the previous result with the naive 
approximation for a Lorentzian density profile 



A0 



I approx 



4M 1 — ro/r arctan(r/ro) „ , ^ 

; 7 — ; — 7 lo^ b < R 

it 1 — ro/R arctan(it/ro) 



(22) 



In Fig. 0, we plot the exact result in Eq. (^) and the naive approximation 
in Eq. (^2]). Here, we have taken M = Mcaiaxy = lO^^M© and R = -Roaiaxy = 
100 kpc as an example for galaxies. In this case, with tq = 0.2i?Gaiaxy, the 
maximal net defiection occurs at 6 ^ Svq and is around 15% bigger than the 
value at 6 = -Rcaiaxy? while for the naive approximation the lensing effect 
inside the galactic halo is always smaller than the effect at 6 = -Rcaiaxy 
Again, we think that a detailed calculation is convenient. 

As a summary of this section, we plot in Fig. ^ the normalized net defiec- 
tion A(j)/A(j){b = i?) as a function of the normalized impact parameter b/R 
for the three specific density profiles considered along the analysis: constant 
density, Gaussian and Lorentzian distribution densities. Such a normaliza- 
tion allows for a clear comparison among the three profiles and is independent 
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Figure 2: The same as in Fig. [^ but for a Lorentzian density profile p{r) = 

j^^^. M = Mcala^y = IO^^Mq, R = i?Gala.y = 100 kpC and Tq = 0.2i?Gala.y 

are used. 

of the mass and physical radius of the compact object. Only for the case of 
stars (Gaussian profile) the lensing effect at the interior of the star is substan- 
tially amplified with respect to the outside effect (the inside effect should be 
compared with the effect aX h = R). For galactic halos, the maximal inside 
net deflection is slightly bigger than aX b = R, while for an object of constant 
density the inside lensing effect is always smaller than at b = R. 

Finally, we plot the focal length / as a function of the impact parameter 
b. The focal length is deflned as the distance at which the lens focuses the 
signal (see Sec. § for details) 



A0 



m • 



(23) 



In Fig. 1^, the (normalized) focal lengths for the three different proflles are 
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Figure 3: Normalized net deflection A(j)/A(f){b = i?) as a function of the 
normalized impact parameter b/R for a constant density profile (solid line), 
a Gaussian density profile (dashed line) and a Lorentzian density profile 
(dotted line). In the last two cases, Tq is taken to be tq = 0.2R. 
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(24) 



We postpone to Sec. § (where focal lengths are discussed in detail) the com- 
ments about the quality of the different gravitational lenses considered here. 
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Figure 4: Normalized focal length f/f{b 
impact parameter b/R. 



R) as a function of the normalized 



3 Neutrino absorption 

While neutrinos are able to travel across the massive object, we have to 
consider their interactions with the matter inside; these will indeed reduce 
the neutrino flux and thus the efficiency of the signal amplification gained 
through lensing. At the energies of interest for cosmic neutrinos {E > 1 — 
lOGeV), we have deep inelastic scattering off individual nucleons. Charged 
currents will cause the conversion of neutrinos into charged leptons, which 
will later be absorbed or decay, while neutral interactions will deflect them 
by angles in general large compared to the lensing effect. In both cases, the 
interacting neutrinos will be lost for lensing ampliflcation. 

Calling (7 the interaction cross-section (which depends on the neutrino 
energy) and N{x) the number density of scatterers along the neutrino path, 
we have the well-known formula for an inflnitesimal flux variation due to 
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interaction with matter: 

d^ = N{x).a{E^).^.dx (25) 

The neutrino flux attenuation after passing through a layer of matter of depth 
X is thus given by: 

^x, E,) = <l>oe-'^(^'') lo ^('-')d-' (26) 

Here we deal with objects having a spherical symmetry, so it is easier to 
express Eq. (^6]) in terms of the radial variable r = yxM-P^, b being the 
impact parameter of the neutrino respect to the center of the object. We 
then get 

-2a(£.)/^^7V(r)- 



<^{b,E^) = %e "'"■''' "V^^^ (27) 

where R is the radius of the object. To solve this equation, we need to collect 
information on the relevant cross-sections, the possible mass density profiles 
and the composition of the object (that is, how to relate mass density to 
number density). 



3.1 Neutrino cross-sections 

At current accelerator energies, neutrino-proton and neutrino-neutron scat- 
tering differ (as the valence quark content is different), as do neutrino and 
antineutrino interactions. From experimental results [0, 0, we adopt the 
following approximations : 



aZ = 0.89 a'^/c ^Z = 1-70 a'Jc 
aZ = 1.70 a^^ aZ = 1-03 a'J'c 



(28) 



using the tables of values for neutrino-isoscalar nucleon from Quigg, Reno et 
al. i 

Above 10^ GeV all these interactions become identical as sea quarks 
largely dominate. We may thus directly use the values for scattering on 
an isoscalar nucleon (again from 0) : 

(^tot = (^tot = (^tot (29) 
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Strictly speaking we should interpolate between those two situations using 
appropriate structure functions, however, in view of the results, the above 
approximations are sufficient. 

Due to the smallness of the ratio — ^^ — , neutrino-electron interactions 
will be neglected compared to neutrino-nucleon ones, except in the particular 
case of the Glashow resonance z/gC" — > W~ —>■ anything, that occurs at an 
energy Ep^ =6.3 PeV. At that energy we have [Q 

Having taken into account the various components of the cross section, we 
rewrite the formula (pTf ) giving the transmission probability of the neutrino 
as: 

fJrl VnT 

PT{h,E^) = exp[-2Y.a^^{E^) / iVx(r) ^== ] (31) 

where X stands for each kind of scatterer. The number densities Nx{r) de- 
pend on the composition and mass density profile of the lens in the following 
obvious way : 

Nx{r) = ^ p{r) (32) 

where A4x and Mx are respectively the fraction in mass and the molar mass 
of constituent X; p is the mass density of the lens. Let us now investigate the 
probability of transmission for the two particular cases of physical relevance 
that have been discussed in the previous section, i.e. stars and galaxies : 

3.2 Neutrino transmission through a star 

All the following calculations apply in a good approximation to all the Sun- 
like stars 0. We will suppose that the star is made up of 25 % of Helium and 
75 % of Hydrogen in number; we then get the following relations between 
number and mass densities : 

NHe{r) = 0.142 mo/ (j^) (33) 

Nnir) = 0.427 mo/ (^) (34) 

Igr 



^Neutrino flux attenuation in the Sun has also been considered in 
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that we used in our numerical calculations. We adopted for the mass den- 
sity profile a Gaussian distribution (normalized to the total mass M, in the 
approximation tq < OAR) that reads: 

I M _(^^2 , . 

in function of the total mass M and radius R of the star, for a distribution 
in which most of the object's mass is concentrated in a sphere of radius rg 
(we have normalized to the total mass). We performed detailed calculations 

R. 

5 ' 

250 M , „„^MHe . „mM, 



in the case tq = -f ; using Eqs. (B^l) and (BB) the transmission probability 



Pt = ■§;; (|3T|) then becomes 

Prib, E^) = exp { — [a"""' — — + a""" ——' 

^ _/5r-,2 rdr 
b \/r'^ — 62 



r 25 M er/(5 Jl - ^) ^, 

exp { — '-— e '^r) 

n R^ erf {5) 



X 



We present in Fig. ^ and ^ the plots of the transmission probability of 
neutrinos across the Sun as a function of h/ Rq] we have added the density 
profile on the same figure. For antineutrinos we have drawn the curve at 
Ey^ = 6.3 PeV where Uee~ is dominant; in this particular case, transmission 
probability can be written as 



25 M erf {5^ 



fe2 



PTib,Q.3PeV) = exp{ — ^^ . . "^ e^^n^ x 

TT R^ erf [5) 

,..-(6.3Pe,.)|2^ + ^]} (37) 

as electron and proton densities are supposed to be equal in stars. 

As we see from Fig. ^ and P, neutrinos with energies of 1 TeV or more are 
completely absorbed or scattered away in the core of the star. Even at 100 
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Figure 5: Transmission probability in function of b, the impact parameter, 
for a neutrino energy E^ = 10 GeV, 10^ GeV, 10^ GeV, 10^ GeV, 10^ GeV, 
10^ GeV (solid lines, from left to right); the Gaussian density profile p{r) = 



e ^ R^ is the dotted hue. 



GeV the interactions are frequent enough to significantly attenuate the flux 
of neutrinos coming out of the star. Unfortunately, the zone where lensing 
could be the most efficient (as the star's core acts as a "real lens") is thus 
ruled out because of flux attenuation for neutrinos of energy 100 GeV or 
more. But we can see that under 100 GeV, interactions of neutrinos inside 
the star have nearly no incidence on the outcoming flux, which is totally 
recovered after lensing by the star. 

3.3 Neutrino transmission through a galaxy 

Neutrinos passing through a galaxy may interact either with its visible matter 
or with the surrounding halo of massive relic neutrinos. These last interac- 
tions become significant only for incoming neutrinos at ultra high energies 
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Figure 6: Transmission probability in function of b, the impact parameter, 
for an antineutrino energy E^ = 10 GeV, 10^ GeV, 10^ GeV, 10^ GeV, 10^ 
GeV, 10^ GeV (solid lines, from left to right). Dashed line on the right is for 
Ep = 6.3 PeV and corresponds to z/gC" scattering. 

(of the order of lO^^eV) ^ and we may neglect them in our present frame of 
work. 

Concerning the visible part of the Galaxy, a rough estimate gives an av- 
erage density of stars of 1 pc~^, which corresponds to a negligible probability 
for the neutrino to encounter a star during its passage through the galaxy, 
even in the worst case if it traverses the whole disk and the bulge. 

We thus conclude that the passage of neutrinos through the galaxy won't 
decrease their flux, and hence do not affect the lensing effect. 

4 Signal enhancement 
4.1 Geometry 
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4.1.1 General relations 

While a significant enhancement will only be reached for well - aligned lens, 
detector and source, we need to consider unaligned elements when dealing 
with extended sources or when taking into account the detector size. We 
follow closely [E|, and remind the basic geometric relations 



(3 = 6 



D 



LS 



D 



-a 



so 



(38) 



where a is the deflection angle; /3, the angular position of the source 
and 6*, the angular position of the image. See Fig. (|^). 
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Figure 7: The geometry of the lensing event. The lens is located at a 
distance Dql from the observer; the source, at a distance Dso- The angular 
separation between the lens and the source is an angle (3 and the position of 
the image is at an angle 9. The deflection angle is a. 

This relation is equivalent to: 



or 



'So far, it has been called A0. 



/5 

(3Dlo 



-f 






(39) 
(40) 
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where ^o, ^o are the values of 9, b in the case of perfect ahgnment. 
The deflection, the characteristic length and angle obey: 





a = 


'-fgib) 




bo = 


OqDlo 




Oo = 


bo 
Dlo 


where 2M ■- 


2GM 
~ c2 ~ 


■ Rschw. 


its mass). 







AMgih) ■ ^^^g^ (41) 



Dso 
Dls 



is the Schwartzschild radius of the lens (M is 

If Rl < b {Rl is the radius of the lens), we are in the usual (photon) case 
(OUTcase): the deflection is given by 4M/b and g{b) = 1. 

As we have seen in section 2, in general g{b) ^ 1 ,and explicit expressions 
were given for Newtonian objects, assuming either Gaussian or Lorentzian 
density profiles. 

An interesting quantity is the focal length / of the lens for given b: 



.m-(^y (42) 

A constant f{b) would signal a perfect lens, where a plane wave focuses in 
one point. Solving graphically / = Dlo determines whether lensing occurs, 
and what its quality is. 

4.1.2 Alignment 

To observe neutrino lensing will usually require huge amplification, and we 
will see that this only happens in the case of perfect alignment {[3 = 0). An 
interesting situation occurs when, in a region of b, say [bmin, bmax], 



g{b) = ^b' (43) 



/4M7 ■ %^ ^ 1 (44) 

Dso 
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In this case indeed the deflection angle is proportional to b, the focal 
length is nearly constant and the object becomes a good lens: all the neutri- 
nos in that region of b are focused on the same point. (The second equation 
simply puts the observer at the focus). 

In this situation, all the particles passing through a ring [bmm, bmax] or 
even the full disk if bmin = focus on the detector. The signal enhancement 
is then given by the ratio of the area of this ring or disk to the area of 
the corresponding disk in case no lensing happens. For a Gaussian density 
profile, the region [0, tq] has a reasonably good lens behavior. Slightly outside 
of the focus area, either a ring or disk would in principle be observed, but 
this requires an angular resolution and statistics unrealistic for neutrinos; our 
only hope is to notice the enhancement factor. [^ 

4.2 Amplification factor. 

4.2.1 General relations. 

The most general relation for the magnification (here, the signal enhance- 
ment) /i is 0: 



_ AQ 



det^. 
39 



A 



OiAOi 



(5A(5 



(45) 



where Afi, Afig stand for the solid angle that the particle trajectories 
span in the sky, respectively with and without lensing effect; 

/?, 6 are the angular coordinates of the source and the image; 

A, Aq are the surfaces of the source and the image referred in a same 
plane (for instance source or lens plane); 

The last equality holds because (3 and 6 are aligned. The index i refers 
to a specific image (in general there is more than one). 

The calculation of the magnification depends on the relation between j3 
and e (given by Eqs. (|8|,|lD). 

^°It is in principle possible to obtain more than one ring, if /(&) crosses / = Dlq niore 
than once, although this does not happen for the Gaussian and the Lorentzian density 
profiles. 



In the OUTcase, we have two images and the resulting magnification is 
well known: 



/^± = 7 y + ^^ ^ ° ±2 (46) 

/^ — ' 1^ = f^OUT (48) 

For the INcase, no general relation holds but we will consider two simpli- 
fied cases (linear and quadratic): 



If g{b) = ^b\ a = AM^b => ^^^ = {^] =^ f^Lk = const. (49) 

A Gaussian density profile, taken close to the center is similar to the 
quadratic case. Note the surprising result that the magnification is indepen- 
dent of j3. As g{b) = '-fb'^ only between 6 = and b = tq (for the Gaussian 
profile), the relation holds for (3 < tq/Dlq- The well known OUTcase be- 
havior is recovered for P ^ tq/ D^o- 



ligib)=Sb,a = 4MS =, ^. = 1 + ^^ = 1 + 1^1 (^0) 

The linear case occurs near b = tq with both Gaussian and Lorentzian 
profiles. There is only a significant magnification near alignment. Otherwise, 
the lens only refracts the signal: it just changes the apparent position of the 
source. 

4.2.2 Alignment. 

Extended sources. This far, we have only considered point-like sources 
and detectors, and as a result, their magnification appears to diverge when 
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source, lens and detector are aligned Eq. (^G]). For an extended circular 
source of uniform brightness and of angular radius Os, the magnification 
factor (OUTcase) is §: 



/^Sc/T = ^ I'' 2n(3d(3^iiP) (51) 

J el + 4^2 

= ^ (52) 



9s 



For the INcase, one evaluates the integral above in the same way, using 
the relations ( ^ , p| ) . For the case g{b) = '-fb'^, i.e. a = AMjb, the image 



of the source (in the lens plane for instance) is a full disk of angular radius 
^disk = ro/Dio- The amplification factor is then |^ 



S -^ ( "disk \ fr-A\ 

For the case a = 4MS, during alignment, the image of the source is a 
ring of radius ^o and width ^ 6s- In the limit 6s ^ 0, the amplification is 
the same as in the OUTcase: 



/i^ = — = 1 + ^ 05^0 260 .. 

Ao 66s 6s 



Extended detector. Finally, we discuss the case of point-like source aligned 
with the lens and the observer. At first glance, the magnification, /x^ or ix^ 
diverges. However, we have not yet taken into account the finite size of the 
detector. The geometry of our problem is symmetric: one can exchange /3 for 



^^ When evaluating the integral, one must substitute /g ^ by Jj^ '*"'=' V^^ _ This takes 
into account the finite size of the lens: the relation a = AM'jb holds only in the interval 

[0',&disk]- 
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^0 (the angle LSO), Dls for Dlo and 9s for O0, the latter being the angular 
radius of the detector seen from the source. The amplification then reads: 



l^OUT 



l^disk 



26^0 Dlo 
6*0 Dls 
I ddisk Dlo 
V O0 Dls. 



(56) 
(57) 



The factor Dlo/Dls is due to the fact that the angle 9^ is not viewed 
from the observer but from the source. When we evaluate the magnification, 
i.e. the surface of the image with and without lensing effect, we have to refer 
to a same plane. Doing so leads to the factor Dlo/Dls- 

Resulting amplification. Until now we have preferred a description with 
angles. Maybe the situation is clearer if we translate the expressions into 
distances in a given reference plane 




Figure 8: To the angular language corresponds formulations in distances 
referred to a same plane. Here, we refer in the lens plane. Ex: The source 
has an angular separation from the lens (3, which corresponds to a distance 
1] in the source plane or equivalently 1]^ in the lens plane. 



Choosing to evaluate quantities in the lens plane, 6*0 — > OqDlo 
SsDlo = -Rs, the radius of the source in the lens plane, /? -^ PDlo 
the position of the source in the lens plane, and OqDls/Dlo 



bo,Os 



e^D 



LS 
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the radius of the detector in the lens plane, as seen from the source. The 
magnification then reads: 



bo_ i,S — 2feo ,,«*_____ — 2bo 

S 



I^OUT — Tf^ f^OUT — 'M f^OUT — R 



2 / \ 2 



,,/3 _ f bo\ i,S _ ( bdisk 1 1,0 _ I bjusk. 

h^disk - U^^ l^disk - y rL J ^J'disk ~ y RL 



(58) 



As we know by Eq. (145|), the magnification factor is given by the ratio 
-^. Which formula to apply results from a competition between the three 
angles P, 9s and O^Dls/Dlo, or, equivalently, between the lengths r]^, Rg 
and R^. Discarding immediately the off- alignment case, we must compare 
two situations, using the lens plane images. If Rg (resp. R^) dominates, 
the image in the lens plane is a ring of width Rg (resp. -R^). Finally, if 
the image is a disk, it is independent of the source and the detector during 
the lensing but the area of the image without lensing is TcRg (resp. tiR^ ). 
So, the answer is: check which length [r]^, Rg or R^) is largest and use the 
corresponding expression for the magnification factor (resp. fi^, fi'^ or fi'^). 
Notice the magnification factor is inversely proportional to the characteristic 
length, or its square, so that the smallest value for the magnification factor 
is always selected. 

Note that in all cases we have the following requirements, in order to use 
the approximated expressions: 



P,es,e, « 9, (59) 

V^R^s.R^ « bo (60) 

4.3 Approximations 

4.3.1 The distant source: Dso ~ Dsl > Dlq 

If the source is far away, the expression for the impact parameter and the 



magnification /i«^ simplifies: bo ~ ^jAMg{h) Dlo, do ^ y^^fff^, R^ ^ R^, 

i/qjjj, ~ ^, /i^.^^ ~ (^) The following inequalities hold: r]^ -C ??, i?^ < Rs- 

A consequence is that in this case one should take care for i?0. Indeed, it 
is possible that R^ -C R0; then the magnification fj,^ applies. 
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4.3.2 Binary systems: Dso ~ Dlq ^ Dls 

If the source and the lens are close to each other, as for binary systems (one 
object being the lens; the other, the source 0), the impact parameter and the 

magnification /i^ simplifies: bo - J4Mg{b)DLs, do - \/^^iP^, Rg ^ R. 



V iO 

2 



^^ouT - Ws ^^d ^^disk - (^) • We have also: t]^ ^ 7] and R^ <^ R^. 

As i?0 < 1 km and as the source is usually macroscopic Rs ^ 1 km, 
it is never needed to take into account the finite size of the detector. So 
ix'^ never applies. For binary systems the most dangerous approximation is 
Rs "^ bo. This should be checked. If this is not realized, one should use the 
full expressions, Eqs. ( ^ , pl| ) . 

5 Applications 

We will know concentrate on possible applications and a few examples. In 
practice, we first test if we are in an OUTcase or not. Then we determine 
the expected shape of the image (discussion of the intersection of /(6) with 
/ = Dlo and its characteristics. We estimate the magnification by deter- 
mining the largest length between r]^, Rg and R^. We check all these lengths 
are smaller than bo to know if we may apply the approximations for the mag- 
nification. Finally, we compute the expected signal for a given detector (in 
our case Super-Kamiokande, SK, or IceCube). 

5.1 The Sun 

We consider the Sun as a first example. 

Neutrinos passing outside the Sun don't focus close enough to Earth to 
be of use; in fact for them to focus on Earth, we would need 

bo - 30.000km ~ 5%Rq (61) 

which is clearly well inside the Sun, so the OUT solution does not provide 
any sizeable effect. 

-"^^Henceforth "binary systems" will point the case Dlo ^ Dls 
-"^^For a treatment of binary systems, see Ref. M 
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We turn to the INcase, using as announced a Gaussian density profile as a 
reasonable approximation of the matter distribution. Unfortunately, the plot 
of the focal length never crosses the line / = D^o = 1 a.u. The smallest value 
is at about 22 a.u., which means Uranians can perform wonderful neutrino 
lensing experiments using the Sun as lens. For them, any source would be 
amplified in turn as the Sun sweeps in front of it! (Jupiter cannot replace 
the Sun as a useful lens for us, as its mass is about IO'^Mq). 

Even if we are not at the focal point, there is yet some amplification due 
to the improved convergence of the beam. It provides: 

The effect is significant if D^o — f\ say D^o ^ [0-3/; 3/], to have a magnifi- 
cation higher than 2. Beyond 3/, this last effect is negligible. 

The Sun is thus not a good candidate amplifying a neutrino signal. Even 
if the neutrinos are seriously defiected, the defiection itself is in practice 
impossible to detect, as it has been done for photons, this is both due to the 
insufficient detection capacity and to the limited angular resolution: 1° in 
the best cases, compared to defiections of less than one arcminute. 

5.2 Stars 

It is well known that stars are good lens candidates. Indeed, the OUTcase 
test tells us: 

6o^ 20i?o(^)^(^)^ if/^Lo«/^L5 (63) 

^o~ 3' — VMoy 

where M is the mass of the lens and Mq, Rq are the solar mass and 
radius. 

For distant sources, since the closest star is yet at about 1.3 pc, the 
OUTcase applies for lensing on stars. The phenomenology is the same as 
for photons. If the source is aligned, the image is a ring of radius 6o- The 
magnification is given by: 
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^Q \Mq) V pc 
30.000km m^(ft^)^ ifDz.o»/^L5 (64) 



,'o„ « I3xl0«(te)(^)'(flf)' iiR,»R-; ^g^j 

^^ovT " 39,000 (f)(JL)^(^)i(^) ifB.«flJ 

The first line applies when the size of the detector determines the radius 
of the ring (this never happens if the source is a star or a galaxy); the second 
is typical for a star or a galaxy as source. However, the "normal" neutrino 
signal of a star (i.e., assumed to be similar to the thermonuclear reactions 
in the Sun) cannot be detected through such amplification. In the most 
favorable case (short D^o and Dso distances): 



N,. 



'events- U ^ ) (, ^^^ j ^'OUT ^e 




/la.u.\^ g 



M y /100pc\ / Ipc 
V Dso ) VA 



0,9X10- I -I (-11-1 l^Jil^l A'. 







The situation could be different if for some reason stellar-size objects hap- 
pened to be strong emitters of energetic neutrinos. In the case of stellar 
lenses, galaxies are in general too wide to be used as a neutrino sources: the 
magnification is trifling if the angular radius of the galaxy 6s is larger than 
the angular radius of the ring Oq 0. 

5.3 Binary Systems 

We discuss know the INcase, which holds for binary systems. We will con- 
centrate on the most promising situation, i.e. a good lens: as we have seen 
before, the image of the source is then a full disk. For our estimations, we 
have chosen a Gaussian density profile with tq = 0.2Rl, the latter giving the 
typical size of the disk. The magnification is thus: 



/^L. = (^) =(i^) (66) 




25 \R 



(67) 



^"^Os is typically larger than 10 ^ for a galaxy and 6*0 is smaller than 10 ^ for a star. 
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The source should thus be smaller than the lens by at least one order of 
magnitude. 

The very question here is to know if the disk is focused on the observer. 
If the focal length at tq is too short, the disk effect will be negligible, and the 
effect will be dominated by an outer part of the lensing star (ring situation); 
if it is too long, no other part of the lensing star will be able to focus on the 
signal on the observer, and the magnification is given by Eq. (|62D . Notice 
here a = AM^yb'^ ^ 4M(6/ro)^ and Dlo ~ Dso- We have the focal length: 




2 

/ Mq\ /la.u. 



\M J \D 



'LS 



(68) 



Clearly, the result f disk/ Dlo ~ 1 is achievable. This situation is really 
promising, providing the source is small enough compared to the lens, so that 
the magnification is significant. Binary stars will however seldom meet all 
the conditions; more exotic systems, with a compact and intense neutrino 
source, are needed. 

5.4 Galaxies 

We now turn to another class of lens candidates, namely galaxies or rather 
the galactic halos. If the source is also galaxy, which is an extended object, 
we will always be in the INcase. Indeed, from Eq. (^TJ), 

. 1 , .1 



6o<V^SI^«14kpo(^)'(§^)' (69) 

Now, the radius of galaxies is typically of this length but, since we will 
focus on the dark matter halo of the galaxy, the INcase applies in most 
situations. 

iFiom. Fig. H, the maximum deflection angle stays close to AM/Rl. More 
precisely, it is always less than 7r/2 ■ AM/Rl. So, 



.„„ < 3 X 10- (JL\ (1^) (70) 



W^^MqI \ Rt 
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and, in order to focus on Earth, 



O,, > /„,. = J±. « 30GPC (^] C-^] (^] (71) 

Since the radius of the Universe is only about 5 Gpc, the neutrinos must 
pass rather near the center of the galaxy, often through its visible part. They 
survive however since the probability to meet a star is tiny, (at the difference 
of photons, which are absorbed by dust clouds). 

For the Lorentzian density profile, the problem is that we have no precise 
idea of the value of ro, except that it is small enough to explain the velocity 
dispersion curves inside galaxies. Here, we discuss only the case 60 > ''"o- The 
deflection angle is then nearly constant: a ~ AM/ Rl ■ [7r/2 — (7r/2 — 1)6]. 
The image is a ring and the magnification is, providing R^ -C -Rf : 



/I00kpc\ / M V' flOMpcY f DsoV' f Di-^^ 



^ "^ ^' V Rs jW'Mq) [ Dlo ) VlGpcj l^lGpcj 

The magnification and the number of events are too tiny to be observed. 
As in the previous section, a small, intense neutrino source is preferable. 



6 Conclusions 

We have studied in some detail the lensing of neutrinos through models of 
astrophysical objects. The main difference with photons rests in the possibil- 
ity for medium-energy neutrinos to cross even the core of stars, and this has 
the important result that the quality of the lens is greatly improved. This 
would have been very promising if the focal length of the center of the Sun 
happened to coincide with the radius of the Earth orbit, unfortunately this is 
far to be the case, and the focusing occurs closer to Uranus. Provided small 
and energetic sources exist, binary systems, where one of the companions 
is the emitter and the other a larger star, or galaxies could provide sizable 
enhancements of the signal. 
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Appendix 

In this appendix, we present the procedure for calculating the deflection of 
neutrinos from straight-line motion as they pass through a gravitational field 
produced by a compact object of mass M. Here, we only consider in detail 
the case concerning the INside solution, i.e. when the neutrino flux passes 
through the object. The procedure is shown for the three specific density 
profiles explained in Sec. |^: constant density, Gaussian and Lorentzian dis- 
tribution densities. 

The procedure starts with the calculation of the mass and pressure of the 
object at a radius r, m{r) and p(r) respectively, for a given density profile 
p{r). We also calculate $(r) and the metric as explained in the main text. 
Second, we derive from the metric the equation of the orbit (for the case of 
constant density we show in detail the derivation) and the smallest radius 
reached by the particle. Finally, we calculate the net defiection. 

We begin by calculating the net defiection for the case of a constant 
density profile p(r) = p 



m{r) = fpr^ = M (^j^y 

A 2M r2 f. 2lvr 



p{r) = p 



n ft TM /-I 2M 7-2 






2V R 2V R R^ 



-\_'2M_( r_Y 
r \r) 

(72) 
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For a massless particle 



Po 



-E 



P dX 

P^ = L 



y — 


/ 




1 


/-I 2M r2 

/!- — 7^ 




3 /, 2M 1 

2V 1 — ir~2\ 




Pr = 


^ 1 dr 

i.m^^ydx 






, P'^ = 


d4, _ L 
dX r2 







(73) 



and the equation of the orbit {p^ = 0) is then 



f, iw 1 L 2I, 
vi — ir-2 yi — IT 



2i\/ 7-2 



R^ 



dr 



1-- 



3J\// 1 
' 2R3 I? 



du 



:^-#)-«^(i-¥ 



/ 1 1(3 /I 2ir 1 /1 2M 7- 



(74) 



Integrating Eq. ( 

0In(m) 



gives 

00 + ^(1 



1 — "la ) + arcsin 



b_ f-i _ M 



+ 



2_R _R2 



1 - ^1 - v/r^F^yi 



i?2 b 



6m 



- -csin{A [1-11(1-31,)]} 



(75) 



+ arcsin 



{6m 



3M 
2/? 



0-*)]} 



where the integration constant is defined so as (pm^r = R) = (t>ovT{r = R) 
(with 00 the initial incoming direction). The particle reaches its smallest r 
0: 



when ^ 

dX 



dr 
dX 



then 











1 


fe2 




1- 


-^c 














[(,. 




1 /-I 2A/ r2 \ 
-2V^ fl li^j 


\ 


1- 


2M 
R 









(76) 



Ur 



1 + 3M f 1 _ _fcL 



+ OiM^u^) 



>IN(n = Mr 



.o + f + ^(l-Vl-f)+ arcsin 



b_ (-1 _ M 



+ liVl^---n{A[l-||(l-3t)]} 



(77) 
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and the net deflection is 



A0 



4M 
b 

4M 



iib> R 



/?2 



2 arcsin 



(i-f) 



I 3M b /I 



ij2 



2arcsin{A[l-||(l-^)]} if 6 < i? 



4M 



if 6>i? 



^ if 6 = i? 







if 6 = 



(78) 
where the outside solution is also included for completeness. 

Next, we present the analysis for the Gaussian and Lorentzian distribution 
densities. As explained in the main text, in both cases, it is possible to make 
some approximations (see Eqs. (|15|,|16D in Sec. H) that allow us to derive 
an analytical expression for the net deflection. For a Gaussian distribution 
density proflle p(r) = poe~^ ^^° 



mlr] 



]^2^ {R^-r^)/rl ro/r e'-''''-O0F/2crf(r/ro)-l 
^ vol Re"'"-'" 



^ '^'■o0F/2crf(R/ro)-l 



M rolre^ ''''o v/i^/2erf(r/ro)-l 
-^ ro/-Re-"^''''o0F/2erf(iJ/ro)-l 

2A/ TQ/re^ ^''O V7r/2erf(r/ro)-l 
ro/R 



R - i"-^^/'''o^/2cTi{R/ro)-l 



(79) 



dt" 



==f, <|)(r) = 
=^ ds^ = - 

V ^ ro/Re'^ /'■o^/2erf{H/ro)-l/ 

where the error function is deflned as erf (2;) = -^ Jq dt e~* . The equation of 
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the orbit is 



dr 



M 



g{fl2-r.2)/.2_^ 



^ ro/Ke-^^/''Ov^/2erf{i?/ro)-l 



^\i^-7^ 



1- 



p2 /^2 



2M 



dy \ ^ R 



0F/2yroerf(l/yro)-e-l/(2"'o)" 



r-o/R e^^ /''O 75^/2 erf (_R/ro )- 1 



i?2/r2 
2AJ rp/re ^/2 erf (»-/ro)-l 

-R _r2/^2 

rg/fle ' 0^/2orf(fl/ro)-l. 



"2„,2-i 



^~y' 



0{M^u 



where y is defined as 



(80) 



M ro/re^'/''o0F/2erf (r/ro) - 1 
" ^ro/i?e^'/'^o0F/2erf (i?/ro) - 1_ 

Integrating Eq. (|^) gives 

0in(z/) = 00 + ^ (1 - Vl-;&) + arcsin(6|/) 



+ 



2M 



ro/Re^^/''o^AF/2 



^ ro/Re-^^/''o0F/2crf(/?/ro)-l 

X 1^1 -;f erf (i?/ro) - Vl-6Verf (1/yr 
e-"/'o erf(^l-f^]-erf 







y/I^ 



1 I ^/n^2„,2\ 



The particle reaches its smallest r at ?/max = i + 0{M^u^), then 



0(l/ = i)= 0o + f + ^(l-;/r^)+^ 



X 



c?2 /^2 



/iie-" /'•oV;F/2 



(81) 



(82) 



r-o/-Re^^''''o0F/2erf(Ryro)-l 



1 - ;^ erf (i?/ro) - e-'''/n^erf ( Jl 



62 K 



/?2 



'"0 



(83) 
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and the net deflection is 



( AM_ 
b 



A0 



^fl-A/l-^) + ^^ 



i?2 



ru/Re" /''Uv^/2 



^ ro/-Re-"^''''o0F/2crf(R/r())-l 



X 

4M 
6 

4M 
R 

4M 



'"0 



1 _ ^ erf (i?/ro) - e-^Vn?erf (^y^l - ;^^ 

if 6>i? 
if 6 = i? 



^fi"/--, 



'Ov/^/2crf(j;/ro) 



527;^ 



ro/Re^ /'■Ov^/2erf(_R/ro)-l'"" 



if 6<i? 
if 6 = 



For a Lorentzian distribution density profile p{r) = j 

^ m{r\ — M '• l~''o/'' arctan(r/r()) 

=^ lliyi ) — IVl ^ 1-ro/R arctan(i?/ro) 

1-ro/r arctan(r/ro)-| log! -5^ 
^ V / i? l-ro/R arctan(iJ/ro) 






( 



ds^ = - 



2M ^^''''■''1''' ^^=*^''(''/''o)"l '°s(ii5^) 
° l—ro/R arctan(_R/ro) 



i? 



+ fl + 2M l-r„/rarctan(r/n,) \ ^^2 ^ ^2^^2 _ 
\ R 1— ro/-H arctan(i(/ro) / 

The equation of the orbit is 



dr 



M 



2 I 2 

1 loo-/ '' +''0 



-R 1—ro/R arctan(_R/ro) 



if 6>i? 



if 6<i? 



rft^ 



(84) 



(85) 



dy 



2M 1—roy arctan(l/roy) 
R l—ro/R arctan(i?/ro) 



/ r2 + r2 \ 
l-ro/r arctan(r/ro)-^ log I „ \^ J 
]^ 2M V" +''0 / 



^-^^^ 



I — tq/R arctan(J?,/rQ) 

0(MV) , 



(86) 
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where y is defined as 



1 

y = - 

r 



jVf 1 - ro/r arctan(r/ro) - \ log (^^^f) 
~ 1 — tq/ R arctan(_R/ro) 



R 



(87) 



Integrating Eq. (^) gives 



?IN 



{y) 



+ ^ 1 - Jl - -^ + arcsin(6|/) 



2M 



1 



.10 



b 1—ro/R arctan(i?/ro) R 



< y 1 — "la arctan(/?/ro) — VI — fe^y^ arctan(l/roy) 



1 + -^ 



arctan 



62 
R^ R 



62 ro 



arctan 
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y^ 1 



62 roy 



The particle reaches its smallest r at ymax = t + O(M^-u^), then 



0(y = i)= 0o+f + ^(i-yr^) 



2M 



1 



.m 



6 1— ro/-R arctan(i?/ro) R 



i?2 



and the net deflection is 
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arctan(i?/ro) 



I 2 

1 + ^ arctan 



62 
R^ R 



1+ 



62 T'O 



A0 



if 6> i? 



1-a/1 



i?2 
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1 + ;3i2 arctan 
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